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CN| ■ Quantum mechanics (QM) disfavours local realism because it violates Bell's theorem [T1I2] . a theorem 
that must be valid for all local hidden variable theories (LHVT). Experimental tests of local realism 
^ c~| in high energy particle physics have been proposed and carried out (3]-[5], but no decisive conclusion 
Q_i has been drawn |6j. We find that the correlated distribution between polarization vectors of two 
■ entangled mesons predicted by QM can be reproduced in a natural way, by interpreting the two-body 
_ decay of vector mesons as a measurement to their polarization vector. The interpretation provides a 
. method to simulate the correlated decay of entangled vector meson pairs, and a potential experimental 
approach for the discrimination between QM and LHVT in particle physics. Related data analysis can 
. be carried out at currently running experiments and could give deeper insight for the understanding 
!> of local realism. 

Debates over the interpretation of QM have lasted for decades. In their famous paper in 1935 (7j , Einstein, Podolsky 
and Rosen (EPR) argued that QM is incomplete by considering an Gedankenexperiment, which is now known as the 
EPR paradox. The paradox challenges the principle of uncertainty that the position and momentum of a particle 
can't be precisely measured simultaneously. This paradox also extends to other pair of conjugate physics quantities. 
In 1951, Bohm expressed the paradox with particle spins in his book 8J and considered a set of alternative theories 
retaining locality and reality, which are called LHVT [9]. It had been considered that one of the LHVT may replace 
QM, but Bell concluded that no LHVT can reproduce all the predictions of QM [TJH], since the predictions of LHVT 
would satisfy the Bell's inequality while those of QM violate it. Other forms of inequalities similar to the original 
Bell's inequality are derived PHE] for the experimental discrimination of the LHVT and QM. 

Since 1972, a series of EPR experiments have been carried out using the entangle photon pairs [12ill6| . and their 
results favored QM. But conclusions can not been drawn decisively because these experiments could not be treated 
loophole free [17]. Tests of the Bell's inequality in particle physics have also been considered. Back to the 1960's, 
the entanglement of K°K° pair was noticed and studied J3J- Violation of a Bell's inequality have been observed in 
the B°B° mesons from T(4S') decay by the Belle collaboration [HE]. But it was argued that the inequalities used is 
not a genuine Bell's inequality so that it can't make a discrimination between QM and the local realism theories [6]. 
Tornqvist suggested that the weakly decay particles, such as the AA pairs generated from the decay of J/ip or r] c , 
could be used for the test of nonlocality of the predictions of QM by assuming the CP invariance [18]. The DM2 
collaboration has studied this with about 10 3 J/ip — > AA events [IS], but the statistics is not sufficient to give a 
decisive conclusion. Test of the local realism with the vector meson pairs from the pseudoscalar cascade decays, such 
as ?y c -> VV -> (PP)(PP), was also suggested [20][21]. 
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In this letter, we try to provide a straightforward understanding of the polarization correlation of entangled vector 
meson pairs. In r] c 's rest frame, the cascade decay process of r\ c to vector mesons V\Vi is illustrated in Fig. [TJ The 
z-axis is chosen to be parallel to meson Vi's momentum. The two vector mesons form an entangled state, and their 
polarization vectors are not fixed. After the decay of meson V\, its polarization is given by the momentum projection 
of its daughter particle d\ on the x-y plane. Another daughter particle of V\ is not illustrated as it has the opposite 
momentum component on the x-y plane to that of d±. The meson V^'s polarization is determined in the same way. 

For the entangled state of two vector mesons, the total spin is S = 1 according to the conservation of the C and 
P parity. In the situation that the total spin projection s on the z-axis is 0, both mesons are transversely polarized, 
and the wave function of the system is [2TjIl2"T] 

|*) = -^(|1),|-1),-|-1),|1),), (1) 

where the subscript z means the z-axis. It can also be expressed in terms of bases in the x-y plane as 

|*) = ^=(|0) Q jO) tt -|0) Q |0) Q J, (2) 

where the subscript a means an arbitrary axis in the x-y plane, and aj_ indicates another axis on the same plane and 
perpendicular to a. 

The correlated distribution between the polarization vectors of the mesons is the same as the probability that the 
spin of one meson has zero component along arbitrary axis a\ and another meson has zero component along arbitrary 
axis OL2- It can be calculated directly with QM. The result is given in Ref. [20] and reads, 



P(m ai =0,m Q2 =0) = L(0U 2 (0|*)| 2 

cx -sm 2 ( ai ~a 2 ) (3) 

= - sin 2 Aq, 
2 

The same form of distribution can be obtained in another way with the understanding of the entanglement of 
polarization. When one meson decays, its polarization vector is determined and has an azimuthal angle a. According 
to Eq. [2] the polarization vector of the second meson should have the azimuthal angle of a ± tt/2. That is similar 
to the case that when we measure the spin of one of the entangled spin- 1/2 particles in one direction, another one 
should have a definite spin value in the same direction. However, the polarization vector of the second meson is finally 
determined by measurement, i.e. its own decay. 

The relation between the polarization vectors and the decays of vector mesons are also illustrated in Fig. [TJ The 
polarization vector Pi of meson V% lays on the projection of the meson's daughter particle's momentum on the x-y 
plane. We denote the axis with the same direction as X, which is the a axis in Eq. [2] The a± axis of Va is denoted 
as Y, which is perpendicular to both the z and X axis. The projection of the momentum vector of V^'s daughter 
particle on the x-y plane is denoted as P2. The angle Aq between X and P2 is exactly the same angle between the 
polarization vectors Pi and P2. As both the vectors are perpendicular to the z-axis, the angle is also equal to the 
angle between the decay planes of the two vector mesons. 

It is obvious that the probability distribution of angle Aa can be calculated if the angular distribution of P2 is 
known. The angular distribution of the final state particles in the decay can be described by the helicity formalism [22] , 
which is widely used in particle physics. For the two-body decay, as the two final state particles are back-to-back in 
their center-of-mass frame, their directions can be characterized by the polar angle 9 and azimuthal angle ip. The 
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cross section of a spin-1 vector meson decay into two spin-0 particles can be written down as 

da (6, tp) da 3 



dfl 



dcos 



47T 



(4) 



where -D M0 



is the Wigner D-Matrix D J mX with spin j = 1, and spin-projection m — M and A = 0, and A is a matrix 
element without any angular dependence. Yy(9, tp) is the spherical harmonic function with 1=1. M depends on the 
selection of spin quantization axis. For the decay of V2, the spin quantization axis is chosen to be in the x-y plane 
with m = 0, which is the Y-axis, so the cross section has an angular dependence of |Y" 1 O (0', tp')\ 2 oc cos 2 6', where the 
polar angle 9' is defined as the angle between the momentum and the Y-axis, and the azimuthal angle tp' is defined in 
the X-z plane. 

The probability distribution should be expressed with 9 and tp in the original xyz coordinate frame. It can be 
obtained with a simple transformation: 



P(cos9,tp) = P(cos6',tp') 



d(cos8',tp') 



d(cos9, tp) 



P(cos9',p'), 



(5) 



where the Jacobian determinant is equal to 1 due to the rotation invariance. Then the probability of Pi and P2 have 
an angle Aa can be obtained by integrating the angular distribution above over 9: 



P(Aa) = J 



P(cos 9, tp)d cos ( 



dp 



dAc 



P(cos9, tp)dcos9, 



(6) 



where Aa and tp differ by a constant, which depends on the selection of the coordinate system. 

It is clear that the value of P(Aa) is irrelevant to the direction of the Y-axis. To simplify the calculation, the 
Y-axis is renamed z'-axis, and chosen to coincide with the x-axis without loss of generality. The y'-axis is the same 
as the y-axis and the a:'-axis is anti-parallel to the z-axis. With this definition, the relations between angles (9, tp) and 
(9' , tp') can be written down: 

cos tp' sin 9' = —cos 9, 

sin tp' sin 9' = sin tp sin 9, (7) 
cos 9' = cos tpsva.9. 



So that the angular distribution is given by 



P(cos 9, tp) oc sin 2 9 cos 2 tp. 



This definition also gives tp = Aa — tt/2, so that the distribution between the polarization vectors is 

P(Aa) oc cos 2 tp = sin 2 Aa. 



(8) 



(9) 



Eq. [S] gives the same form of result as that given in Eq. |3 The prediction on the correlated distribution of QM is 
reproduced naturally. 

The calculation above provides a method to simulate the quantum entanglement in the process of cascade decay 
of T] c —> VV —> (PP)(PP). The key of the simulation is to sample the directions of final state particles. That is 
straightforward from Fig.[T] At first, either of the two vector mesons is selected and denoted as V\. In the rest frame of 
rjc, the decay of V\ would follow the angular distribution of |Y 1 ±1 (6', tp)\ 2 oc sin 2 9. The angles 9 and tp can be sampled 
from this distribution. The sampled p determines the quantization axis of the second vector meson V2. The decay of 
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V-i follows the angular distribution of \Yi(6', (p')\ 2 oc cos 2 9', from which the angles 8' and <p' are sampled. By rotating 
the unit vector constructed with the sampled 8' and (p' back into the same coordinate frame of V\ , the decay angle 8 
and ip of V2 can be obtained. 

This method is implemented for the process of ?/ c —»•(/></>—>• K + K~ K + K~ . To check whether the simulation 
method works as expected, 10,000 events are generated, and the correlated distribution f(Aa) is extracted. The 
result is plotted in Fig. [5J The distribution is fitted with the function of p ■ sin 2 A<x It is clear that the distribution 
satisfies that given in Eq. [3] and Eq. [9j The method works well. 

Another more important conclusion of the understanding is that it provides a potential approach to discriminate 
QM and LHVT in high energy experiments. An analogy between the decay and the polarimeter used in the optical 
experiments to test the Bell inequalities could be made. Though the direction of the "polarimeter" can't be selected 
at will, high statistical data make the revelation of the correlated distribution of the polarization vectors possible. 
On the contrary, it is difficult to obtain such correlated distributions in the optical experiments, because one need 
to place the two polarimctcrs in all possible direction and try all the combination of directions. A variation of the 
CHSH inequality has been derived in Refs. [2TJII2T] . The inequality is violated by QM at some points due to the special 
correlated distribution predicted by QM. In other words, if the predicted correlated distribution is observed, one could 
say definitely that the inequality would be violated, and QM would be favored. 

It is possible to measure such a correlated distribution in current high energy experiments. The branching fraction 
of J/ip -t jr) c is about 1.7 ±0.4%, of n c -> <j><p is about 2.7 ±0.9 x 10~ 3 and cj> -> K+K~ is about 48.9 ±0.5% [23]. At 
the BES-III experiments, lO 9 J/-0 particles are expected to be recorded in the next few years [23]. About 10 5 events 
for the rj c — > <fxj> — > K + K~ K + K~ process are expected to be collected. The statistics would be sufficient for the 
measurement of the correlated distribution. The analysis result might be of fundamental importance. 

In summary, with the understanding of the polarization correlation, we reproduce the same form of correlated 
distribution predicted by QM of the polarization vectors of entangled vector meson pairs. This provides a way for the 
simulation of the correlation in the pseudoscalar cascade decay. The measured correlated distribution in high energy 
particle experiments can be used to make a discrimination between QM and LHVT. The measurement is possible in 
existing particle physics experiments. 
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Figure 1: Decay of a pair of entangled vector mesons from the decay of rj c . Pi and P2 are the polarization vector of 
meson V\ and V2, respectively. 
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Figure 2: The correlated distribution of the polarization vectors of the two vector mesons. The x-axis is the angle 
between the polarization vectors. The distribution follows the formula of f(Aa) oc sin 2 Aa. 
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